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Abstract 

Motivated by the study of conditional stability of traveling waves, 
we give an elementary center stable manifold construction for 
quasilinear parabolic PDE, sidestepping apparently delicate regular- 
ity issues by the combination of a carefully chosen implicit fixed-point 
scheme and straightforward time-weighted energy estimates. As an 
application, we show conditional stability of Lax- or undercompressive 
shock waves of general quasilinear parabolic systems of conservation 
laws by a pointwise stability analysis on the center stable manifold. 
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1 Introduction 

In this paper, extending our previous work in the semilinear case [Z5], we 
show by an elementary argument that an asymptotically constant stationary 
solution 

(1.1) u{x,t)=u{x), l'u(x) - < Ce-^l^l, 

^ > 0, of a general quasilinear second-order parabolic system 

(1.2) ut = J^{u) := h{u)uxx - h{u, u^), 

X, t e M, li, /i e M", b e M"^", possesses a local center stable manifold with 
respect to H^. Combining this result with ideas from [ZH, HZ, RZ, TZ2, 
TZ3, Z5] we then establish conditional stability of nonclassical viscous shock 
solutions of strictly parabolic systems of conservation laws, similarly as was 
done in [Z5] for Lax shocks of semilinear parabolic systems. 
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1.1 Existence of center stable manifolds 

We briefly describe our results. FoUowing [Z5], assume: 

(AO) h e C"=+i, k>2. 
(Al) ^a{b) >e>0. 

(A2) The linearized operator L = about u has p unstable (positive 

real part) eigenvalues, with the rest of its spectrum of nonpositive real part. 
(A3) \diu{x)\ < Ce^^l^l, ^ > 0, for 1 < j < A; + 2. 

Then, we have the following basic version of the Center Stable Manifold 
Theorem. 

Proposition 1.1. Under assumptions (A0)-(A2), there exists in an 
neighborhood of u a Lipschitz (with respect to H"^) center stable manifold 
M-cs, tangent to quadratic order at u to the center stable subspace T^cs of L 
in the sense that 



for u e M-cs where lies o^nd 11^, denote the center-stable and stable eigenpro- 
jections of L, that is (locally) invariant under the forward time- evolution of 
(1.2) and contains all solutions that remain sufficiently close to u in forward 
time. In general it is not unique. 

Combining with ideas of [TZl, Z5], we readily obtain by the same tech- 
nique the following improved version respecting the underlying translation- 
invariance of (1.2), a property that is important in the applications [Z5]. 

Theorem 1.2. Under assumptions (A0)-(A2), there exists in an neigh- 
borhood of the set of translates of u a translation invariant Lipschitz ( with 
respect to H"^) center stable manifold M.cs, tangent to quadratic order at u to 
the center stable subspace of L in the sense of (1.3), that is (locally) in- 
variant under the forward time- evolution of (1.2) and contains all solutions 
that remain bounded and sufficiently close to a translate ofu in forward time. 
In general it is not unique. 

1.2 Conditional stability of nonclassical shocks 

Now, speciahze to the case 

(1.4) h{u, u^) = f{u)^ - {db{u)u^)u^ 




(1.3) 
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that (1.2) corresponds to a parabolic system of conservation laws in standard 
form 

(1.5) ut + f{u)^ = {h{u)u^)^, 

u, f eW,be M"''", X , i e R. Following [ZH, HZ], assume: 
(HO) /, b G C"=+2, k>2. 
(HI) ^a{b) >9>0. 

(H2) A± := df{u±) have simple, real, nonzero eigenvalues. 

(H3) ^a{dfii - b\i\^){u±) < -^ICP, ^ > 0, for all ^ e R. 

(H4) Nearby u, the set of all solutions (1.1) connecting the same values 
u± forms a smooth manifold a &U C M^, m° = m. 

(H5) The dimensions of the unstable subspace of df{uJ) and the stable 
subspace of df{u+) sum to either n + £, i = 1 (pure Lax case), n + £, £ > 1 
(pure overcompressive case), or < n with £ — 1 (pure undercompressive case). 

Assume further the following spectral genericity conditions. 

(Dl) L has no nonzero imaginary eigenvalues. 

(D2) The Evans function D{X) associated with L vanishes at A = to 
precisely order £. 

Here, the Evans function, as defined in [ZH, GZ] denotes a certain Wron- 
skian associated with eigenvalue ODE {L — \)w = 0, whose zeros correspond 
in location and multiplicity with the eigenvalues of L. For history and ba- 
sic properties of the Evans function, see, e.g. [AGJ, PW, GZ, MaZl] and 
references therein. 

As discussed in [ZH, MaZl], (D2) corresponds in the absence of a spectral 
gap to a generalized notion of simplicity of the embedded eigenvalue A = 
of L. Thus, (D1)-(D2) together correspond to the assumption that there 
are no additional (usual or generalized) eigenvalues on the imaginary axis 
other than the transational eigenvalue at A = 0; that is, the shock is not in 
transition between different degrees of stability, but has stability properties 
that are insensitive to small variations in parameters. 

With these assumptions, we obtain our remaining results characterizing 
the stability properties of u. Denoting by 

(1.6) af < < ■ ■ ■ < 

the eigenvalues of the limiting convection matrices A± := df{u±), define 

(1.7) e(x, t) := ^ (1 + i)-V2e-|a:-a-t|VMt + ^ (1 + ^)-l/2g-k-a+t|VMt^ 

aj<0 at>0 
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a7<0 



+ x(x, V (1 + kl + t)-'/'(l + \x- an\)-'/\ 



2^ 

at>0 



and 
(1.9) 



^2(a;, t) := (1 - x(a;, t))(l + |a; - art| + ti/2)-3/2 

+ (l-x(x,^))(l + k-a+^|+^^/')-'/^ 



where x{x^ t) = 1 for x G [a^ t, a+t] and zero otherwise, and M > is a 
sufficiently large constant. 

Proposition 1.3. Conditions (H0)-(H5) imply (A0)-(A3), so that there ex- 
ists a translation-invariant center stable manifold M.cs ofu and its translates. 

Theorem 1.4. Under (H0)-(H5) and (D1)-(D2), u is nonlinearly phase- 
asymptotically orhitally stable under sufficiently small perturbations vq G 
lying on the codimension p center stable manifold Aics ofu and its translates 
with 1(1 + \x\'^Y^'^vq{x)\h'^ < Eq sufficiently small, where p is the number of 
unstable eigenvalues of L, in the sense that, for some a{-), a*, 

\dl{u{x,t)-u^'+^^'\x)) \ < CEo{e + i^i + i^2){x,t), 

\u{;t)-u'"-+''^'^\H^<CEo{l+t)--^, 

(1-10) < CEo, 

\a{t)\<CEo{l + t)-'/\ 
\a(t)\<CEo(l + t)-\ 

< r < 1; where u denotes the solution of (1.5) with initial data Uq — u-\-Vq 
and u'^ is as in (H4). Moreover, u is orbitally unstable with respect to small 

perturbations not lying in M., in the sense that the corresponding solution 
leaves a fixed-radius neighborhood of the set of translates of u in finite time. 

Remark 1.5. Pointwise bound (1.10) yields as a corollary the sharp L*' 
decay rate 

(1.11) \u{x,t)-u'^''+<'^\x)\Lv <CEo{l + t)-'^^^-p\ l<p<oo 

obtained by the LP (rather than pointwise) analysis of [Z5]. However, we 
obtain here the additional information that phase a* + a approaches a limit 
time-asymptotically at rate t~^. 
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Remark 1.6. Mixed over- or undercompressive type shocks are also possible 
[ZH, HZ], though we do not know of any physical examples. These can be 
treated with further effort as described in [HZ], Sections 5-6, or [RZ]. 

1.3 Discussion and open problems 

As discussed in [GJLS], for semilinear problems like those considered in 
[GJLS, Z5], for which the nonlinear part of the associated evolution equations 
consists of a relatively compact perturbation of the linear part, construction 
of invariant manifolds reduces essentially to verification of a spectral map- 
ping theorem for the linearized flow, after which the construction follows 
by already-well-developed theory. However, for quasilinear equations, the 
usual fixed-point construction of the standard theory docs not close due to 
apparent loss of regularity. This appears to be a general difficulty in the con- 
struction of invariant manifolds for quasilinear systems; see [Li, LPSl, LPS2] 
for further discussion. 

We overcome this in the present case in a very simple way, by (i) intro- 
ducing a carefully chosen implicit fixed-point scheme for which the infinite- 
dimensional center stable part satisfies a standard Cauchy problem forced by 
the finite-dimensional (hence harmless in terms of regularity) unstable part, 
and (ii) making use of a straightforward time-weighted energy estimate for 
the Cauchy problem to control higher derivatives by lower ones, the latter of 
which may be estimated in standard linear fashion. For related arguments, 
see for example [MaZ2, MaZ3, Z2, Z3, RZ, TZ3]. 

Existence of invariant manifolds of quasilinear parabolic systems has been 
treated by quite different methods in [LPSl, LPS2] via a detailed study of 
the smoothing properties of the linearized fiow. The advantage of the present 
method, besides its simplicity, is that carries over in straightforward fashion 
to the compressible Navier-Stokes equations of gas dynamics and MHD [Z6] 
for which the linearized flow, being of hyperbolic-parabolic type, is only 
partially smoothing. A disadvantage is that we do not see by this technique 
how to obtain smoothness of the center stable manifold, but only Lipschitz 
continuity; this seems to be the price of our simple energy-based approach. 

As we are mainly interested in stability, it is not smoothness but quadratic- 
order tangency (1.3) that is our main concern. An observation of possibly 
general use is that this weaker property is satisfled whenever the underlying 
flow is at u, whereas global regularity of the center stable manifold 
would require global C^"*"" regularity, o; > 0, of the flow. Moreover, it is easily 
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verified in the course of tlie standard construction of a Lipschitz invariant 
manifold, and so we obtain this key property with essentially no extra effort. 

Regarding conditional stability, the main novelty of the present analysis 
is that we carry out a pointwise iteration scheme in order to treat nonclas- 
sical shock waves (for further discussion regarding the need for pointwise 
estimates, see [HZ]). It seems an observation (though elementary) of possi- 
bly wider use that our center stable manifold construction can be used to 
obtain pointwise control on the solution in this way, extending a bit further 
the basic ideas of [Z5] . As our study of conditional stability was partly moti- 
vated by discussions in [GZ, AMPZ, Z7] of certain unstable undercompressive 
shocks and their effect on solution structure through metastable behavior, it 
seems desirable to fit such nonclassical waves in the theory. 

An interesting open problem is to investigate conditional stability of a 
planar standing shock u{x,t) = u{xi) of a multidimensional system of con- 
servation laws 



which likewise (by the multidimensional arguments of [Zl, Z2, Z3]) reduces 
to construction of a center stable manifold, in this case involving an infinite- 
dimensional unstable subspace corresponding to essential spectra of the lin- 
earized operator L about the wave. 

Plan of the paper. In Section 2, we establish existence of center stable 
manifolds for general quasihnear parabolic PDE. In Section 3, we estabhsh 
conditional stability on the center stable manifold by a modification of the 
pointwise arguments of [HZ, RZ] in the stable (p = 0) case. 



Defining the perturbation variable v := u — u, we obtain after a brief com- 
putation the nonlinear perturbation equations 
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2 Existence of Center Stable Manifold 



(2.1) 



Vt — Lv — N{v), 



where 



(2.2) 



'X 
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denotes the linearized operator about the wave and 

N{v) := (b{u + v){u + v)^j: - h{u)u^:c - b{u)vj:j: - {db{u)v)% 



(2.3) 

- \ h{u + v,u^ + v^) - h{u, u^) - hu{u, Ux)v - h^^ {u, u^)v^j 

is a quadratic order residual. We seek to construct a C^{H'^) local center 
stable manifold about the equilibrium v = 0, that is, a locally invariant C*^ 
manifold tangent (in Frechet sense) with respect to to the center stable 
subspace Ecs- 

2.1 Preliminary estimates 

Denote by S„ and the unstable and center stable subspaces of L and 
and lies the associated spectral projections. 

Proposition 2.1 ([TZl]). Under assumptions (A0)-(A3), L generates an 
analytic semigroup e^^ satisfying 

(2 d) -^1. - uj , 

\e ii„|L2_^2 < Cu,e ^ , 

for some /3 > 0, and for all cu > 0, for all t > 0. 

Proof. Standard semigroup estimates for second-order elliptic operators; see 
[TZl, Z5], or Appendix A. □ 

(1 Ixl < 1, 
\x\ > 2, 



N\v) :=p(^)Ar(^). 



Lemma 2.2 ([TZl]). Assuming (A0)^(A3), the map : ^ %s C''+^ 
and its Lipschitz norm with respect to v is 0{S) as 5 — > 0. Moreover, 

(2.5) \N\v)\h2 < C\v\l2. 

Proof. See Appendix A. □ 
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Corollary 2.3. Under assumptions (A0)-(A3), 



(2.6) 



for some /? > 0, and for all a; > 0, /or all t>Q, with Lipschitz bounds 



(2.7) 



2.2 Fixed-point iteration scheme 

Applying projections 11^, j — cs,u to the truncated equation 

(2.8) vt- Lv = N\v), 

we obtain using the variation of constants formula equations 

Ujv{t) = e^^'-''-^%v{toj) + I e^^'-''^IijN\v{s)) ds, 

J to,j 

j = cs,u, so long as the solution v exists, with toj arbitrary. Assuming 
growth of at most |t>(t)|jj2 < Ce^* in positive time, we find for j — u 
using bounds (2.4) (ii) and (2.7) (ii) that, as io,M +oo, the first term 

e^*-*~*'''"''n„f (to.u) converges to zero while the second, integral term converges 
It^°° e^^*~*'*nuA^'^(v(s)) ds, so that, denoting w :— UcsV, z :— U^v, we have 

/+00 
e^^'-'^UuN\{w + z){s)) ds. 

Likewise, choosing to,cs = 0, we have 

(2.10) w{t) = e^^UM + f e^^^-'^U^sN\(w + z)(s)) ds, 

Jo 

Wq := Ilcsf (0). On the other hand, we find from the original differential equa- 
tion projected onto the center stable component, after some rearrangement, 
that w satisfies the Cauchy problem 



wt - b{u)wxx = Ilcsb{'"')^xx - n„6(M)m 



XX 



(2-11) 

+ U^sM{w + z) + U^sN^ {w + z) 
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with initial data Wq — llcsv{0) given at t — 0, where 

(2.12) M{v) := {db{u)v)ux - hu{u,Ux)v - hu^{u,Ux)vx. 

We shall use these two representations together to obtain optimal estimates, 
the first for decay, through standard linear semigroup estimates, and the 
second for regularity, through the nonlinear damping estimate (2.5) below. 

Viewing (2.10), or alternatively (2.11), as determining w — yV{z,wo) as 
a function oi z, we seek 2; as a solution of the fixed-point equation 

(2.13) z = f{z,Wo):=r{z,Wiz,Wo)). 

As compared to the standard ODE construction of, e.g., [B, VI, TZl, Z5], in 
which (2.9)-(2.10) together are considered as a fixed-point equation for the 
joint variable {w, z), this amounts to treating w implicitly. This is a standard 
device in situations of limited regularity; see, e.g., [CP, GMWZ, RZ]. 

It remains to show, first, that W, hence T, is well-defined on a space of 
slowly-exponentially-growing functions and, second, that T is contractive on 
that space, determining a solution z = z{wq) similarly as in the usual 
ODE construction. We carry out these steps in the following subsections. 

2.3 Nonlinear energy estimates 

Lemma 2.4 ([Z5]). Under assumptions (A0)-(A3), for all 1 < p < 00, 

< r < 4, 

(2.14) I n„ I i;,p_>vt^r,p , \Tlcs\wr,p-^W^'P ^ C . 

Proof. See Appendix A □ 

Proposition 2.5. Under assumptions (A0)-(A3), for 5 sufficiently small, if 
the solution of (2.11) exists ont & [O;^]? then, for some constants 6*1,2 > 0, 
and allO<t<T, 

(2.15) \wit)\l, < Ce-^^>o|^2 + C f e-''^'-'\\w\l, + \z\l,){s)ds. 

Jo 

Proof. Let us first consider the simpler case that b is uniformly elliptic, 

sR6:=^(6 + 6*) >^>0. 
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Taking the inner product in x of Yi^j=G ^V'^ against (2.11), integrating by 
parts and rearranging the resulting terms, we obtain 

dt\w\l,{t) < -(9>,6(s)a» + C(|«;||. + 1^11,4), 

<-0\dlw\h + C[\w\l:. + \z\l,), 

> 0, for C > sufficiently large, where we have repeatedly used the bounds 
of Lemma 2.2 and of Lemma 2.4 with p = 2 to absorb the error terms coming 
from the righthand side of (2.11), and have used Moser's inequality to bound 
\dxN^ {v)\l2 < C\v\loo\v\h'a < C5\v\h3 whenever does not vanish, so that 
|'y|i/2 < 26. Using the Sobolev interpolation 

\w\l. < C-'\diw\l. + C\w\l. 

for C > sufficiently large, and observing that 

(2.16) < C\z\l^ 

by equivalence of norms on finite-dimensional spaces, we obtain 

dt\w\l,{t) < -e\w\]j, + C {\w\l, + \z\l,) , 

from which (2.15) follows by Gronwall's inequality. 

To treat the general case, we note that 3?(t(6) > ^ > by Lyapunov's 
Lemma implies that there exists a positive definite matrix P such that 
3?(P6) > y, whence, by a partition of unity argument, there exists a smooth 
positive definite matrix-valued function P{u) such that K(P6) > |. Taking 

the inner product in x of Yl'^=o^xP{'^)^x''^ against (2.11), integrating by 
parts and rearranging the resulting terms, we obtain 

dtS{w) < -{dlw,Pb{u)^^w)+C{\w\j,, + \z\j,,), 

(2.17) e 

^ > 0, where S{w) := Ylj=o{'^x'^j is equivalent to | • |^2- By 

Sobolev interpolation and (2.16), we have therefore 

dtS{w) < -9E + C{\w\l, + \z\l,), 
9 > 0, from which (2.15) follows again by Gronwall's inequality. □ 
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Proposition 2.6. Under assumptions (A0)-(A3), if solutions Wi,Zi and 
W2,Z2 of (2.11) exist for allt > 0, then, ford sufficiently small, any constant 
9>0 and some C = 0(9), for all t > Q, 



(2.18) 



/ e '^'{wi - W2\h2{s) ds < C\wo,i - Wo^2\h^ 
Jo 

+ C e-^'{\wi-W2\l2 + \zi-Z2\l2){s)ds, 
Jo 



+00 



(2.19) 



e^^'-^Vl - ^2\H2is) ds < C\wi - W2\Hi{t) 

/ + 00 
e^(*-)(|^, - W2\l2 + \zr - Z2\l.){s) ds 

+ Ce^'^\\\wi-W2\\l2 +\\zi-Z2\\l2 ). 



Proof. Subtracting the equations for wi,zi and W2-iZ2, we obtain, denoting 
w :— wi — W2, z :— z\ — Z2, the equation 

(2 201 ^* ~ K^)wxx = Ilcsb{u)z^cc - Ilub{u)wr,^ 

+ U,sM{w + z) + U,,{N\wi + zi)-N\w2 + Z2)) 

with initial data Wq = Wo,i — Wo,2 at t = 0. 

Performing an version of the energy estimate in the proof of Propo- 
sition 2.5- that is, taking the inner product in x of '^]=odiP{u)di'w 
against (2.20), P as in the proof of Proposition 2.5, integrating by parts, and 
rearranging the resulting terms- we obtain 

(2.21) dt£{w) < ~\dlw\l2 + C{\w\l, + \z\l2), 

a > 0, where £{w) :— '^]=Q{diw, Pb{u)diw) is equivalent to Here, 
we have used in a key way the H'^ — > Lipschitz bound on A'"'^ to bound 

\N\wi + Zi) - N\w2 + Z2)\l2 < CS{\w\h2 + \z\h2), S « 1, 

and thus {\d^w,Pb{u)d^w),N^{wi + Zi) - N^{w2 + Z2))\ < C(|w|^2 + |i||2). 
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By Sobolev interpolation and (2.16), we have therefore 

(2.22) dtS{w) + \w\],2 < C{\w\l2 + \z\l2), 

whence (2.18) follows by Gronwall's inequality, S >0, and |£'(wo)| < C\wo\'jji. 

The first line of (2.19) follows similarly, using \£{w{t))\ < C\w\'jji{t). 
Noting that we could in place of (2.22) have rearranged (2.21) as 

d,£(w)<--£(w) + C{\w\l2 + \z\l2), 

and argued as in Proposition 2.5 to obtain 

\w{t)\]j, < Ce-'''\wo\]j^ + C f e-''^'-'\\w\l2 + \z\l2){s) ds 

Jo 

<Ce-'''\wo\jji+Ce''^\\\w\\L2_ +\\z\\l2_ ), 



(2.23) 



substituting in the first line of (2.19), and estimating similarly 

/+00 
e"-'-'\\w\l. + \z\l2){s)ds < Ce'^'\\H\Ll^ + UhlJ, 

we obtain the second line of (2.19), completing the proof. □ 

Remark 2.7. The absence of a uniform Lipschitz bound on 

prevents us from obtaining a pointwise energy estimate on w like the one 
obtained on w in Proposition 2.5. 

2.4 Basic existence result 

Define now the negatively-weighted sup norm 

||/||_, := supe-^*|/(t)|^,., 
t>o 

noting that |/(t)|j^2 < e^*||/||_g for alH > 0, and denote by the Banach 
space of functions bounded in || • ||_^ norm. Define also the auxiliary norm 

^2 .= supe-''*|/W|L^- 

^ t>0 
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Lemma 2.8. Under assumptions (A0)-(A3), for 2>uj <r] < (5 and 5 > and 
Wq G sufficiently small, for each z G there exists a unique solution 

w W{z,wo) e B-ri of (2.10), (2.11), with 

(2.24) ||«;||_^<C(|«;o|^^2 + 5||z||_^) 
and 

(2.25) ||>V(zi, «;o,i) - W(z2, «;o,2) lU^.^ < C5\\z^ - Z2\\-r, + C||«;o,i - «;o,2||-^. 

Proof, (i) (H^ hound) By short-time if ^-existence theory, and the earher- 
observed fact (2.16), an solution of (2.11) exists and remains bounded 
up to some time T > provided that wq is bounded in if^, whereupon 
(2.15) holds. 

Using now the integral representation (2.10), and applying (2.4) (i), (2.7) (i) 
^ \N'{v(t))\ < 5\v{t)l and \w{t)\H2 < e^l*' ||w||_^,t supo<,<T e-''^|w(s)|,j2, 
we obtain for < t < T that 

\w{t)\L2 < Ce^\'\\wcs\L^ + CS{\\w\\.rj,T + |k||-^) / e'^l^-^'e'^l^l ds, 

Jo 

hence, using a; ± > that |w(t)|L2 < Ce'''*' ^|wcs|l2 + 5(||w||_^^r + ||2;||_^) j . 
Applying (2.15), we then obtain 

\w(t)\H^ < Ce^l*l (\w,s\h^ + 5(\\w\\^r,,T + Ikll-,)) 

for all < t < T, and thus ||w||-r;,T < c(^\wcs\h^ + S{\\w\\_rj^T + W^W-vi^^ 
whence, for S sufficiently small, < C{\wcs\h^ + \\z\\-ri)- Since this 

bound is independent of T, we obtain by continuation global existence of w 
and, letting T — >• oo, (2.24) as claimed. 

(a) (H'^ — > Lipschitz bounds) Now consider a pair of data zi,wq^i and 
Z2,wo^2, and compare the resulting solutions, denoting {z,w,wo) '■— {zi — 
Z2, W1 — W2, wo,i — wo^2)- Using the integral representation (2.10), and applying 
(2.4) (i), (2.7) (i), and the definition of || • ||_^, we obtain for alH > 

(2.26) \w{t)\L2 < Ce'^l*l|wo|L2 + C5 / e^l*-^l(|wjj^2 + \z\H^){s)ds. 

Jo 
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Estimating as before 

CS f e'^^'-'^\z\m{s)ds <C5\\z\\-r, f e^^'-'^e'^' ds < CS\\z\\^r,e''' 
Jo Jo 

and, by (2.18) with 6 = Slj together with the Cauchy-Schwarz inequaUty, 
C6 e"(*-^) \w\h^ (s) ds<C6(^ g-'^^*"^)) ( e^^^*"^) \w\l, (s) ds^ 

we obtain, substituting in (2.26), 
(2.27) 

<e^''*(cKl^i + C^5^(||^|li._^ + Pir-.)). 
This yields 

||w(t)||i2^ < C\wo\-r^ + C6{\\w\\l2_^ + ||i||-r,), 

from which (2.25) follows by smallness of 6. □ 

Proof of Proposition 1.1. (i) (Boundedness on a ball) Again, working in 
recall that 

/ + 0O 
e^(*-^)n„7V'5((^^ + ^)(s))ds, 

where w = W{zai)q). Using (2.9), and applying (2.4)(ii), (2.7)(ii), \N\v{t))\ < 
5\v{t)\, and \z(t)\H^ < e"!*! ||^||_^ and, by (2.24), 

\w{t)\H^ < e^l*l||«;||_, < Ce^\'\{\wo\H^+5\\z\\_,), 

we obtain for <t that 

/+00 

hence, using f3±r] > and taking S and juiolj/z sufficiently small, that T(-, wq) 
maps the ball B{0,r) to itself, for r > arbitrarily small but fixed. 
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(ii) (Lipschitz bounds) Similarly, we find using (2.25), (2.7)(ii), and Lemma 
2.2 that 

\\'T'(zi,Wo,l) -f(z2,Wo,2)\\-r, 

/ + 00 



< supCi5M|i||_^ + e-''* / ef^^^-'^\w\H^{s)ds]. 

^ *J t 



Using the Cauchy-Schwarz inequahty and (2.19) with 9 — P to estimate 

-oo , r+oo ^ I /.;, , r+oo 



1/2 



/+0O . r+co , 1/2 / r+oo 

< Cs{\w\Ut) + 1^ e^('-^\\w\l. + \z\h){s)ds) 

< Cae^'M |wo|hi + ||w||l2 + ||i||L2 ) 
and applying (2.25), we obtain 

(2.29) \\f{zi,Wo,i)-f{z2,Wo,2)\\-v < CS{\\wo\\-r, + \\z\\LiJ. 

This yields at once contractivity on 5(0, r), hence existence of a unique fixed 
point z = Z{wo), and Lipshitz continuity of Z from to by the 
Banach Fixed- Point Theorem with Lipschitz dependence on parameter wq. 
(ii) (Existence of a Lipschitz invariant manifold) Defining 

r+oo 

(2.30) ^wo) Z{wo)\t=o = - / e-'^m^N^vis)) ds, 

Jo 

we obtain a Lipschitz function from S^, whose graph over 5(0, r) 

is the invariant manifold of solutions of (2.8) growing at exponential rate 
\'^{t)\ < C'C* in forward time. Prom the latter characterization, we obtain 
evidently invariance in forward time. Since the truncated equations (2.8) 
agree with the original PDE so long as solutions remain small in H^, this 
gives local invariance with respect to (1.1) as well. By uniqueness of fixed 
point solutions, we have Z{0) = and thus $cs(0) = 0, so that the invariant 
manifold passes through the origin. Likewise, any bounded, sufficiently small 
solution of (1.1) in is a bounded small soluton of (2.8) as well, so by 
uniqueness is contained in the center stable manifold. 
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'cs)\h'^ 







(a) (Quadratic-order tangency) By (2.30), (2.4), (2.5), and (2.24), 

e-^m^N^w + z)(s) ds 

r+oo 

<C e-^\\W\]j. + \Z\jj,){s)ds 

< Ci{\m% + < C,{M'h. + \\Z\\%). 

By Z{0) = and Lipshitz continuity of Z, we have < C\wcs\h^, 

whence 

(2.31) mWcs)\H^<Cs\w,s\ji2, 

verifying quadratic-order tangency at the origin. □ 
2.5 Translation-invariance 

We conclude by indicating briefly how to recover translation-invariance of the 
center stable manifold, following [Z5, TZl]. Differentiating with respect to 
a the relation == J-'{u{x + a)), we recover the standard fact that (j) '■— 
is an zero eigenfunction of L, by the assumed decay of Ux- 
Define orthogonal projections 

(2.32) u^-.^^r^, ni:=id-n2, 

onto the range of right zero-eigenfunction := {d/dx)u of L and its orthog- 
onal complement (j)^ in L^, where (•, •) denotes standard inner product. 

Lemma 2.9. Under the assumed regularity h e C'''^^, k > 2, Uj, j — 1,2 
are bounded as operators from to itself for Q < s <k -\-2. 

Proof. Immediate, by the assumed decay oi (f) — and derivatives. □ 

Proof of Theorem 1.2. Introducing the shifted perturbation variable 

(2.33) v{x, t) :— u{x + a{t),t) — u{x) 
we obtain the modified nonlinear perturbation equation 

(2.34) dtv^ Lv + N{v)-dta{(l) + dxv), 
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where L :— ^{u) and A?^ as in (3.17) is a quadratic-order Taylor remainder. 

Choosing dta so as to cancel 112 of the righthand side of (2.34), we obtain 
the reduced equations 

(2.35) dtv 
and 

(2.36) dta 

V G (f)-^, := {(p,v)\(f)\jj2, of the same regularity as the original equations. 
Here, we have imphcitly chosen a{0) so that v{0) e 0-*-, or 

{(f), uq{x + a) — u{x)) = 0. 

Assuming that Uo lies in a sufficiently small tube about the set of translates 
of M, or Uq{x) = {u + w){x — with \w\h2 sufficiently small, this can be done 
in a unique way such that a := a — P is small, as determined implicitly by 

Q = G{w, a) := (0, u{- + a) - u{-)) + (0, w(- + a)), 

an application of the Implicit Function Theorem noting that 

d^gio,o)^{<f>,dM-))^\<l>\l^^o. 

With this choice, translation invariance under our construction is clear, with 
translation corresponding to a constant shift in a that is preserved for all 
time. 

Clearly, (2.36) is well-defined so long as \dxv\L^ < C\v\h2 remains small, 
hence we may solve the v equation independently of a, determining a- 
behavior afterward to determine the full solution 

u{x, t) — u{x — cii{t)) + v{x — a{t),t). 

Moreover, it is easily seen (by the block-triangular structure of L with re- 
spect to this decomposition) that the linear part ITiL = niLIli of the v- 
equation possesses all spectrum of L apart from the zero eigenvalue asso- 
ciated with eigenfunction 0. Thus, we have effectively projected out this 
zero-eigenfunction, and with it the group symmetry of translation. 

We may therefore construct the center stable manifold for the reduced 
equation (2.35), automatically obtaining translation-invariance when we ex- 
tend to the full evolution using (2.36). See [TZl, Z5] for further details. □ 



^Ui{Lv + N{v)) 
_ n^jLv + N{v)) 
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2.6 Application to viscous shock waves 

Proof of Proposition 1.3. Clearly, (HO) implies (AO) by the form (1.4) of /i, 
while (HI) and (Al) are identical. Plugging u — u{x) into (1.5), we obtain 
the standing-wave ODE f{u)x — ib{'^)'^x)x, or, integrating from — oo to x, 
the first-order system 

b{u)u^ = fiu) - f{u_). 

Linearizing about the assumed critical points u± yields linearized systems 
Wt = df{u±)w, from which we see that u± are nondegenerate rest points by 
(H2), as a consequence of which (A3) follows by standard ODE theory [Co]. 

Finally, hnearizing PDE (1.2) about the constant solutions u = u±, we 
obtain Wf = Lj^w := —df{u±)wx — b{u±)w,,.j.. By Fourier transform, the 
limiting operators L± have spectra \f{k) = a{—ikA^ — b±k'^), where the 
Fourier wave-number k runs over all of M; in particular, L± have spectra 
of nonpositive real part by (H3). By a standard result of Henry [He], the 
essential spectrum of L lies to the left of the rightmost boundary of the 
spectra of L±, hence we may conclude that the essential spectrum of L is 
entirely nonpositive. As the spectra of L to the right of the essential spectrum 
by sectoriality of L, consists of finitely many discrete eigenvalues, this means 
that the spectra of L with positive real part consists of p unstable eigenvalues, 
for some p, verifying (A2). □ 

3 Conditional stability analysis 

From now on we specialize to the case (1.4), (1.5) of a system of viscous 
conservation laws Ut = J^{u) = {b{u)ux)x — f{u)x-i with associated linearized 
operator L — dTu{u). 

3.1 Projector bounds 

Let H„ denote the eigenprojection of L onto its unstable subspacc S„, and 
Hcs = Id — H„ the eigenprojection onto its center stable subspace 'Lcs- 

Lemma 3.1 ([Z5]). Assuming (HO)-(Hl), 

(3.1) n^dx = dxiij 
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for j — u, cs and, for all 1 < p < oo, < r < 4, 

\^cs\wr,p^W'''P , \^cs\wr,p^W^-P ^ C, 

|(l + |a;|2)in„/|H4<|/|i., 
|(l + |xnin„/|^4<|(l + |a;ni/|^4. 

|neJ(a;)|<Ce-''l^lsup|e-''l^l/(l/)|. 
y 

Proof. Recalling that L has at most finitely many unstable eigenvalues, we 
find that Uu may be expressed as 

where (pj, j = l,...p are generalized right eigenfunctions of L associated 
with unstable eigenvalues A^, satisfying the generalized eigenvalue equation 
(L — XjY^(l)j — 0, rj > 1, and are generahzed left eigenfunctions. Noting 
that L is divergence form, and that \j ^ 0, wc may integrate {L — XjY^(f)j = 
over M to obtain A^-^ J (j)jdx = and thus f (pjdx = 0. Noting that (t)j, (pj and 
derivatives decay exponentially by standard theory [He, ZH, MaZl], we find 
that 

(f)j = d^^j 

with $j and derivatives exponentially decaying, hence 

3 

Estimating 

3 j 

for 1/p + l/q = 1 and similarly for d^Huf, wc obtain the claimed bounds on 
n^, and Uu, from which the bounds on lies = Id — n„ and lies = Id — n„ 
follow immediately Bounds (3.3) and (3.4) follow similarly. □ 



(3.2) 

and 
(3.3) 

Moreover, 
(3.4) 
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3.2 Linear estimates 

Let Gcs{x, t; y) :— Ilcse-^^^y{x) denote the Green kernel of the hnearized solu- 
tion operator on the center stable subspace Tics- Then, we have the following 
detailed pointwise bounds established in [TZ2, MaZl]. 

Proposition 3.2 ([TZ2, MaZl]). Assuming (H0)-(H5), (D1)-(D2), the cen- 
ter stable Green function may he decomposed as Gcs — E + G, where 

(3.5) E{x,t;y) = J2^P^ 

j=l j \a=a* 

h (^) - -Mtw^)) 

for y <0 and symmetrically for y >0, with 

(3.7) \1%\<C, \{d/dy)l%\<C^e-^\y\, 
and 

\d'^G{x,t;y)\ < Ce-''(l^-^l+*) 

n 

+ ^(r t + e-^i^i) ( J2 t-h-^^-y-'^^'^'^^'e-"^^ 

k=l 

(3.8) + Y: X(|„-,>M>t-^/^e-(-'^7(*-l^/'^^l))V-*e--^ 

>0, aj <0 

I V V .-1/2 -{x-a+it-\y/a- my Mt -nx-\ 

a~>Q,a+>Q 

(3.9) 

\d'JyG{x,t;y)\ < Ce-''(l^-^l+*) 

n 

+ Ct-'^t-i + e-^l^l + 76-^1^1) ( Yl rV^e-(^-^-«^"*)'/^*e-''^+ 

k=l 

+ 2^ ^{|»^*I>|2/|}^ e ^ e 

+ Ml«fct|>l2/|}^ e ^ y 

aj^>0,at>0 
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for < s < 1, for y < and symmetrically for y > 0, for some rj, C, 
M > 0, where are as in (1.6), (5^ > 0, denotes the positive/negative 
part of X, indicator function X{\a~t\>\y\} ^-^ ^ Z^'" Wk^\ — Ivl '^^'^ ^ otherwise, 
and 7 = 1 m the mixed or undercompressive case and in the pure Lax or 
overcompressive case. 

Proof. This follows from the observation [TZ2] that the contour integral (in- 
verse Laplace Transform) representation of Gcs is exactly that for the full 
Green kernel in the stable case p — 0, and that the resolvent kernel satisfies 
the same bounds. Thus, we obtain the stated bounds by the same argu- 
ment used in [MaZl] to bound the full Green kernel in the stable case. See 
Appendix A for further discussion. □ 

Corollary 3.3 ([HZ]). Assuming (H0)-(H5), (D1)-(D2), 
(3.10) 

\ej{y,t) — ej{y, +oo)\ < Cerr/n( ^^^ ), some a > 
X errjn I — = — errjn ' 



y4/3-(t+i); 

\dytejiy,t)\ < C(t-^ + 7t-V2e-''lfl) ^ e-l^+"^"*l'/^* 



for y < 0, and symmetrically for y > 0, where 7 as above is one for under- 
compressive profiles and zero otherwise. 

Proof. Direct computation using definition (3.6); see [MaZl]. □ 
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3.3 Convolution bounds 

Prom the above pointwise bounds, there follow by direct computation the fol- 
lowing convolution estimates established in [HZ] , here stated without proof. 

Lemma 3.4 (Linear estimates [HZ]). Assuming (H0)-(H5), (D1)-(D2), 
(3.11) 

f+OO 



oo 

+00 



oo 

■ oo 



/+0O 
\G{x,t-y)\{l + \y\)-"^dy < C{e + + i^2){x,t), 
-oo 

\G,{x,t-y)\{l + \y\)-^/'dy < C{t--^ + 1)(^ + + ^2)(a;, t), 

/ + 0O 
\e,{y,t)\{l + \y\)-'/' dy <C{l + tr'l\ 
-oo 

/+00 
|e(y,i)|(l + |y|)-^/^dy<C, 
-oo 

/+0O 
|e(y, t) - e(y, +oo) | (1 + \y\)-^'^ dy < C(l + ^)-V^ 
-oo 

for < i < +00; C > 0; where G and e are defined as in Proposition 3.2. 

Lemma 3.5 (Nonlinear estimates [HZ]). Under (H0)-(H5), (D1)-(D2), 
(3.12) 

ft C + OQ 

\Gy{x, t - s; y)\^{y, s) dyds < C{d + + ij2){x, t), 



'0 -'-oo 

"t-i r+oo 



rt—i r+oo 

/ / \Gy^{x,t-s;y)\^{y,s)dyds<C{0 + iPi + iP2){x,t), 

Jo J-oo 

ft r+oo 

/ / \G^{x,t-s;y)\{9 + iPi + iP2){y,s)dyds<C{0 + iPi + iP2){x,t), 

Jt-l J-oo 

pt r+oo 

/ / \eyt{y,t-s)\^{y,s)dyds<C{l + t)-\ 

Jo J-oo 



p+oo r+oo 

/ / \ey{y, +oo)\^{y, s) dy < Cj{l + t)-^l\ 

Jt J-oo 

rt r+oo 

/ / \ey{y,t-s)-ey{y,+oo)\^{y,s)dyds<C{l + t)-^/'', 

Jo J-oo 
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for < i < +00, C > 0, where G and e are as in Proposition 3.2 and 

+ (l + s)-^(^ + ^l+^2)(|/,s). 

We have by standard short-time theory the following additional bound. 
Lemma 3.6 (Commutator bound). Assuming (H0)-(H5), (D1)-(D2), 

/+00 
\{G, + Gy){x,t-y)\{l + \y\) dy < C{0 + i^, + 4^2) {x, t) , 
-00 

for < t < 1, C > 0, where G is as in Proposition 3.2. 

Proof. By standard short-time parametrix bounds [Fr] on the entire Green 

function, we have for < i < 1 that \Gx + Gy\ < Ct~^e mt- is of the same 
order order as G, whence 

/+00 
\{Gx + Gy){x,t;y)\{l + \y\)-'l^dy<C{l + \x\)-''\ 
■00 

By direct computation, integration against l^ajGul, \dyGu\, \dxE\, or \dyE\ 
gives a contribution that is also bounded by C(l -|- yielding the 

result. □ 



3.4 Reduced equations II 

We now restrict to the pure Lax or undercompressive case i = 1, following 
the simple stability argument of Section 3, [HZ]. The pure Lax or over- 
compressive case may be carried out following the similar but slightly more 
complicated argument of Section 4, [HZ]. 

Define similarly as in Section 2.5 the perturbation variable 

(3.15) v{x,t) := u{x + a{t),t) — u{x) 

for u a solution of (1.5), where a is to be specified later in a way appropriate 
for the task at hand. Subtracting the equations for u{x + a{t),t) and u{x), 
we obtain the nonlinear perturbation equation 

(3.16) vt- Lv^ N{v)x - dta{(j) + d^v), 



24 



where L as in (2.2) denotes the hnearized operator about u, (f) — u^, and 

(3.17) N{v) := -(/(« + v)- f{u) - df{u)v) 

where, so long as \v\h'^ (hence \v\l°° and |m|l«>) remains bounded, 

( . Niv) = Oi\v\\v,\), 

^ ■ ^ d,N{v)^0{\v\\v,,\ + \v,\'). 

RecaUing that dxU is a stationary solution of the linearized equations 
Ut — Lu, so that Ldx^O, or 

/oo 
G{x,t;y)ux{y)dy = e^^u^^x) = d^u^x), 
■oo 

we have, applying Duhamel's principle to (3.16), 

/•oc 
G{x,t;y)vo{y) dy 
-oo 

f't f'OO 

~ / Gy{x,t — s;y){N{v) + av){y, s) dy ds + a{t)dxu{x). 

Jo J-oo 

Defining 

/oo 
e{y,t)vQ{y) dy 

+ / ey{y,t- s){N{v) + av){y,s)dyds, 

Jo J-oo 

following [ZH, Z4, MaZ2, MaZ3], where e is defined as in (3.6), and recalling 
the decomposition G = E + 0^ + of (A. 2), we obtain the reduced equations 

/oo 
{Gu + G){x,t;y)voiy) dy 

{Gu + G)y{x,t- s; y) {N{v) + av) {y, s)dy ds, 



-'-oo 

and, differentiating (3.19) with respect to t, and observing that ey{y,s) —r 
as s — > 0, as the difference of approaching heat kernels, 

/oo 
et{y,t)vo{y) dy 

{'^■■^^) rt r+oo 

+ eyt{y,t-s){N{v) + dv){y,s)dyds. 

Jo J-oo 
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We emphasize that this (nonlocal in time) choice of a and the resulting 

reduced equations are different from those of Section 2.5. As discussed further 
in [Go, Z4, MaZ2, MaZ3, Z2], a may be considered in the present context as 
defining a notion of approximate shock location. 

3.5 Nonlinear damping estimate 

Proposition 3.7 ([MaZ3]). Assuming (H0)-(H5), let G H*, and suppose 
that for < t < T, the H"^ norm of v remains bounded by a sufficiently 
small constant, for v as in (3.15) and u a solution of (1.5). Then, for some 
constants 9i^2 > 0, for all <t <T, 

(3.22) \v{t)\l, < Ce-'''\v{0)\l, + C f e-''^'-'\\v\l2 + \a\^){s) ds. 

Jo 

Proof. Energy estimates essentially identical to those in the proof of Propo- 
sition 2.5, but involving the new forcing term a{t)dxu{x). Observing that 
dl{dxu){x) = 0{e~^^^^) is bounded in norm for j < 4, we obtain (in the 
simpler, uniformly elliptic case), taking account of this new contribution, the 
inequality 

dt\v\Ut) < -0\dlv\l. + C{\v\l, + \ait)\') , 

^ > 0, for C > sufficiently large, so long as |w|//4 remains sufficiently small, 
yielding (3.22) as before by Sobolev interpolation followed by Gronwall's 
inequahty. See also [MaZ3, RZ]. □ 

3.6 Short time existence theory 

Lemma 3.8 ([RZ]). Assuming (H0)-(H5), let Mq := |(1 + |x(2)^/So(x)|jj4 < 
oo, and suppose that, for < t < T , the supremum of \a\, and the norm 
of V, determined by (3.16), each remain bounded by some constant C > 0. 
Then there exists some M = M{C) > such that, for all < t < T , 

(3.23) \(l + \x\^f^v(x,t)\l,<Me^'{Mo+ f\aWr)dr). 

Jo 

Proof. This follows by standard Friedrichs symmetrizer estimates carried out 
in the weighted if^ norm. Specifically, making the coordinate change 

v^il + lxiy/^'w, 
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we obtain from (3.16) a similar equation plus lower-order commutator terms, 
and similarly in the equations for for j = 1, . . . , 4. Performing the same 
energy estimates as carried out on (3.16) in the proof of Lemma 3.7, we read- 
ily obtain the result by Gronwall's inequality. We refer to [RZ, Lemma 5.2] 
for further details in the general partially parabolic case. □ 

Remark 3.9. Using Sobolev embeddings and equation (3.16), we see that 
Lemma 3.8 immediately implies that, if \{\ + \x\^Y^^Vq{x)\h4. < oo and if 
k (•) ^)|if4; \a.{t)\ are uniformly hounded onO <t <T, then 

\{l + \x\y/%{x,t)\ and \{l + \x\^)^/\{x,t)\ 

are uniformly bounded on <t <T as well. 

3.7 Proof of nonlinear stability 

Decompose now the nonlinear perturbation v as 

(3.24) v{x,t) ^w{x,t) + z{x,t), 
where 

(3.25) w :— HcsV, z UuV. 

Applying lies to (3.20) and recalling commutator relation (3.1), we obtain 
an equation 



/oo 
G{x,t;y)wo{y) dy 



oo 

t /'OD 

/ Gy{x,t - s;y)Ilcs{N{v) + av){y,s)dyds 

J-oo 

for the flow along the center stable manifold, parametrized by -u; G lies- 

Lemma 3.10. Assuming (HO)-(Hl), forv lying initially on the center stable 
manifold M.cs, 

(3.27) |9^z(x,i)| < Ce-^l^l|«;|^2 

for all < r < 4, some C > 0, so long as \w\h2 remains sufficiently small. 
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Proof. By quadratic-order tangency at u of the center stable manifold to 
Scs, estimate (1.3), we have immediately \z\fj2 < C|?«|^2, whence (3.27) 
follows by equivalence of norms for finite-dimensional vector spaces, applied 
to the p-dimensional subspace E„. (Alternatively, we may see this by direct 
computation using the explicit description of UuV afforded by Lemma 3.1.) 

□ 

Proof of Theorem 1.4- Recalling by Theorem 1.2 that solutions remaining for 
all time in a sufficiently small radius neighborhood jV of the set of translates 
of u lie in the center stable manifold Aics- we obtain trivially that solutions 
not originating in Aics niust exit AT in finite time, verifying the final assertion 
of orbital instability with respect to perturbations not in A4cs- 

Consider now a solution v e A4cs, or, equivalent ly, a solution w e of 
(3.26) with z = ^cs{w) G E„. Define 

C(t):= sup f{\w\ + \w,\){e + iPi+ij2)-\y,t) 

(3.28) 

+ |w|j/4(l + s)3 + |d(s)|(l + s)y 

We shall establish: 

Claim. For alH > for which a solution exists with ( uniformly bounded 
by some fixed, sufficiently small constant, there holds 

(3.29) C(^)<^^2(^o + C(^)')- 



From this result, provided Eq < 1/4C|, we have that ({t) < 2C2E0 
implies ({t) < 2C2E0, and so we may conclude by continuous induction that 

(3.30) C{t) < 2C2£;o 

for all t >0, from which we readily obtain the stated bounds. (By standard 
short-time existence theory together with Remark 3.9, v e if^ exists 
and ( remains continuous so long as ( remains bounded by some uniform 
constant, hence (3.30) is an open condition.) 

Proof of Claim. By (3.3), \{1 + \x\^)^wo\h4 = \{l + \xf)^UcsVo\H* < CEq. 
Likewise, by Lemma 3.10, (3.28), (3.18), and (3.4), for < s < 

.33^^ \ncs{N{v) + av)iy,s)\<Caty^iy,s)l 

^ ■ ^ \{N{v) + av){y,s)\<Cat)'^{y,s)\, 
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^ as in (3.13), while 

(3.32) \U,sdy{N{v) + av){y, s)\ < CC{tfC{l + s)-^e + + ilJ2){x,t). 

Combining (3.31) and (3.32) with representations (3.26) and (3.21) and 
applying Lemmas 3.4 and 3.5, we obtain 



\w{x,t)\< \G{x,t;y)\\wo{y)\dy 

J —oo 

+ \Gy{x,t - s;y)\\Ilcs{N{v) + av){y,s)\dyds 

Jo J-oo 

<Eo r \G{x,t;y)\{l + \y\)-'/'dy 

J —OO 

pt poo 

+ CC{tf / \Gy{x,t-s-y)\^{y,s)dyds 
<C{Eo + atT){0 + ipi + ip2){x,t) 

and, similarly, 

/oo 
\et{y,t)\\uo{y)\dy 
-oo 

+ \eyt{y,t - s)\\flcs{N{v) + av){y,s)\dyds 

Jo J-oo 

^^■^^^ < r Eo\etiy,mi + \y\)-'/'dy 

J —oo 

t P + OD 



+ 



/ / Cat)yyt{y,t-s)\-^{y,s)dyds 

Jo J-oo 



<c{E,+cm{i+t)-\ 
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Likewise, we may estimate for i > 1 
(3.35) 

/oo 
\G^{x,t;y)\\wo{y)\dy 
-oo 



Pt—l POO 

+ / \G^y{x,t-s;y)\\Ucs{N{v) + av){y,s)\dyds 

Jo J -oo 

Pt POO 

+ / \G^{x,t- s;y)\\Ucsdy{N{v) + av){y,s)\dyds 

Jt-l J-oo 

/oo 
\G{x,t;y)\{l + \y\)-'/'dy 
-OO 

p — 1 poo 

+ CC{tf / \G^y{x,t-s;y)\^{y,s)dyds 

Jo J-oo 

pt POO 

+ CCitf / / \G^,{x,t-s;y)\{9 + ^i + ^2) {y, s)dy ds 
J t—i J —00 

<c{Eo + atf){e + A + '4^2){x,t), 

and for < i < 1, substitute for the first term on the righthand side instead 

/ 

J — c 



\{G^ + Gy){x,t;y)\\wo{y)\dy+ / \G{x,t;y)\\dyWo{y)\dy 

J —00 

to obtain the same bound with the aid of (3.14). 
By Lemma 3.10, 

(3.36) \z\H^{t) < C\w\]j2{t) < CC{tf. 

In particular, |;i;|/,2(t) < CC,(tY{l + t)~^ . Applying Proposition 3.7 and using 
(3.33) and (3.34), we thus obtain 

(3.37) \w\H^{t) < C{Eo + C{tf){l + t)-i. 

Combining (3.33), (3.34), (3.35), and (3.37), we obtain (3.29) as claimed. 

As discussed earlier, from (3.29), we obtain by continuous induction 
(3.30), or C < 2C2|wo|Lin//4) whereupon the claimed bounds on \v\, \vx\ and 
\v\Hi follow by (3.33), (3.35) and (3.37) together with (3.27), and the bounds 
on \a\ by (3.34). It remains to establish the bound on \a\, expressing con- 
vergence of phase a to a limiting value q;(-|-oo). 
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By Lemmas 3.4-3.5 together with the previously obtained bounds (3.31) 
and ( < CEq, and the definition (3.28) of the formal limit 



a(+oo) := / e{y,+oo)uo{y)dy 

J —oo 

P + OO P + OO 

+ / ey{y, +oo) {N {v) + av){y,s)dyds 

Jo J-oo 

/oo 
Eo\e{y,+oo)\{l + \y\)-'/Uy 
-oo 

+ / CEo\ey{y, +oo)\^ {y,s)dyds 

Jo J-oo 



< CE, 







is well-defined, as the sum of absolutely convergent integrals. 
Applying Lemmas 3.4-3.5 a final time, we obtain 

/oo 
Hy,t) -e{y, +oo)\\vo{y)\ dy 
-oo 

pt /" + 0O 

+ / / Kiy,t-s) - ey{y, +oo)\ 

Jo J-oo 

X I {N{v) + av) {y, s) \ dy ds 

p+oo P+OO 

+ / \(iy{y,+<>o)\\{Q{v) + av){y,s)\dyds 

Jt J-oo 

/oo 
Eo\e{y,t)-e{y,+oo)\{l+\y\)-^/^dy 
-oo 

pt P + OO 

+ / / \ey{y,t - s) - ey{y,+oo)\CEo'if{y,s)dyds 

Jo J-oo 

P + OO P + OO 

+ / \ey{y,+(x)\CEo'i{y,s)dyds 

J t J-oo 

<CEo{l + t)-'/^ 
establishing the remaining bound and completing the proof. 



□ 
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for pointing out the reference in [GJLS]. 
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A Proofs of miscellaneous lemmas 



We include for completeness the proofs of earlier cited lemmas that were not 
proved in the main body of the text. 

Proof of Proposition 2.1. By sectoriality of L, we have the inverse Laplace 
transform representations 



(A.1) 



where Fes denotes a sectorial contour bounding the center and stable spec- 
trum to the right [Pa], which by (A2) may be taken so that KF^ < a;, 
and r„ denotes a closed curve enclosing the unstable spectrum of L, with 
> /9 > 0. Estimating \e^^^j\ui < |e-^*| |(A-L)~^|i2|dA| using sectorial 

resolvent estimates, we obtain the stated bounds from ^ L^; sec [He]. 
Applying the resolvent formula L(A — L)~^ — A(A — L)~^ — Id, we obtain 



Le^^Uj := [ e^\\-L)-^d\ 



from which we obtain |Le^*n„|i2 < Ce ^* for t < by boundedness of F„, 
yielding the stated bound from L'^ ^ H'^. □ 

Proof of Lemma 2.4. Recalling that L has at most finitely many unstable 
eigenvalues, we find that may be expressed as 

where 4>j, j = l,..-p are generalized right eigenfunctions of L associated 
with unstable eigenvalues Xj, satisfying the generalized eigenvalue equation 
(L — \jY^(f)j = 0, Tj > 1, and (f)j are generalized left eigenfunctions. Noting 
that 0j, (f)j and derivatives decay exponentially by standard theory [He, ZH, 
MaZl], and estimating 
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for 1/p + 1/g = 1, we obtain the claimed bounds on n„, from which the 
bounds on = Id — n„ follow immediately. □ 

Proof of Lemma 2.2. The norm in H'^ is a quadratic form, hence the map 
is smooth, and is as regular as N. The Lipschitz bound follows by 

< 3|^;i-^;2|^^2f sup + sup \dN{v)\LA, 



so that sup|^|^2<5 = 0(5^), sup,^,^2<5 |rfiV(t;)|L2 = 0{6). Finally, 

\N\v)\l2 < |A^(f)|L2 < C\v\]j2 for \v\h2 < 26 by Moser's inequality, while 
N\v) = for \v\h2 > 25, yielding (2.5). □ 

Proof of Proposition 3.2. As observed in [TZ2], it is equivalent to establish 
decomposition 

(A.2) G^Gu + E + G 

for the full Green function G{x,t;y) := e^*dy{x), where 

p 

i=i 

for some constant matrix M G C'*'^^ denotes the Green kernel of the linearized 
solution operator on E„, (pj and (pj right and left generalized eigenfunctions 
associated with unstable eigenvalues Xj, j — 1, . . . ,p. 

The problem of describing the full Green function has been treated in 
[ZH, MaZ3] , starting with the Inverse Laplace Transform representation 

(A.3) G{x, t; y) = e^^5y{x) = j e^\\ - L{e))-^5y{x)d\ , 
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where 

r:= d{X : m < rji - rj2\QX\} 

is an appropriate sectorial contour, r/i, 7^2 > 0; estimating the resolvent kernel 
G\{x,y) := (A — L{e))~^6y{x) using Taylor expansion in A, asymptotic ODE 
techniques in x, y, and judicious decomposition into various scattering, ex- 
cited, and residual modes; then, finally, estimating the contribution of various 
modes to (A. 3) by Riemann saddlepoint (Stationary Phase) method, moving 
contour F to a optimal, "minimax" positions for each mode, depending on 
the values of {x,y,t). 

In the present case, wc may first move F to a contour F' enclosing (to the 
left) all spectra of L except for the p unstable eigenvalues Xj, j = 1, . . . ,p, 
to obtain 

G{x,t;y)^ i e^*(A - L)-^ciA + VResidueA.(e)(e^*(A - L)-^5j/(x)), 

where ResidueA^(£) (e'^*(A — L)^^8y[x)) — Gu{x,t;y), then estimate the re- 
mammg term e^\X- Ly^dX on minimax contours as just described. See 
the proof of Proposition 7.1, [MaZ3], for a detailed discussion of minimax es- 
timates E+G and of Proposition 7.7, [MaZ3] for a complementary discussion 
of residues incurred at eigenvalues in {3?A > 0} \ {0}. See also [TZ2]. □ 
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